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Two recent experiments successfully observed Landau levels in the tunneling spectra of the topo- 
logical insulator Bi2Se3. To mimic the influence of a scanning tunneling microscope tip on the 
Landau levels we solve the two-dimensional Dirac equation in the presence of a localized electro- 
static potential. We find that the STM tip not only shifts the Landau levels, but also suppresses 
for a realistic choice of parameters the negative branch of Landau levels. 

PACS numbers: 73.20.-r, 71.70.Di 
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Topological insulators have recently attracted consid- 
erable interest, especially after their experimental discov- 
ery in two- and three dimensions 1 9 . While insulating in 
the bulk, such materials possess gapless edge states whose 
existence depends on time reversal invarianc c 10 '' 16 l This 
makes the latter robust against time-reversal symmet- 
ric perturbations such as impurity scattering and at the 
same time very sensitive to time-reversal breaking ones 
such as magnetic fields. When the topological insulator 
is a three-dimensional system, the gapless excitations are 
confined to its surface and form a two-dimensional con- 
ductor. 

The Landau quantization of the surface states of the 
three-dimensional topological insulator Bi2Se3 has been 
reported in two recent experiment j 17 l 18 4. In both cases 
the Landau levels were detected using a scanning tunnel- 
ing microscope (STM). The differential tunneling con- 
ductance 31/ dV measures the local density of states of 
electrons at energy eV. In the absence of magnetic fields 
the measured spectra are consistent with a Dirac disper- 
sion of the surface states. In the presence of a magnetic 
field Landau levels are expected at energies 



E n = E 



D 



sgn{n)vp^J 2eh\n\B 



(1) 



where Er> is the energy of the Dirac point, n — 
0, ±1, ±2, ... is the Landau level index, vp is the veloc- 
ity, e is the unit charge, H is the Planck constant, and B 
is the magnetic field. A series of unequally spaced Lan- 
dau levels has been observed in the above-mentioned ex- 
periments, including a £>-independent level at the Dirac 
point. However, while the theory predicts a positive 
(n > 0) and a negative (n < 0) branch of Landau levels 
experimentally only t he po sitive branch has been seen. 
It has been speculate d 17 1 18 1 that the absence of Landau 
levels below the Dirac point may result form an overlap- 
ping of the surface states with the bulk valence band. 
On the other hand, as pointed out in Ref. [IB], near the 
Fermi energy the bulk conduction band overlaps with the 
surface state, but still Landau levels are clearly observed. 
In this paper we investigate another possible reason 
for the absence of the negative branch of Landau lev- 
els, namely the electrostatic effect due to the STM tip: 
The authors of Ref. [17] have noticed that the Dirac 
point in their STM measurements is about 200 meV be- 
low the Fermi level, while the Fermi level determined by 



angular-resolved photoemission spectroscopy is only 120 
meV above the Dirac pointPS, and suggested that this 
discrepancy might be due to the electrostatic interaction 
between the STM tip and the sample. We will elaborate 
this idea further and will demonstrate that such a poten- 
tial may indeed strongly suppress the negative branch of 
the Landau levels. 

In the following we will present the model under in- 
vestigation, sketch the numerical methods and finally we 
will present the results. Close to the Dirac point the sur- 
face states of a topological insulator with a sing le Dirac 
cone can be described by the HamiltoniarPESEI] 



H = v F (p + eA) • a- + V(x) 



(2) 



where p is the two-dimensional momentum operator, A 
is the vector potential and cr are the Pauli matrices. Note 
that cr is in general not the spin operator. For example, 
for Bi2Se3, symmetry requires for the spin operator the 
relation Soce z x <rP2 V(x) is the electrostatic potential 
caused by the tip, which we characterize by its depth and 
width. The depth can be extracted from the experiment 
from the position of the Dirac point. The width is not 
directly known but can be estimated to be of the order 
10 - 20 nmP The results presented later are obtained 
with a Gaussian potential, V(x) = Vgexp(— |x| 2 /2cr 2 ). 

For the numerical treatment it is convenient to intro- 
duce dimensionless units. In the following we measure 
distances in units of an arbitrary length scale xq. The 
energy is measured in units of E$ = hvp/xo, and the 
magnetic field in units of Bq = h/(exo). Numerically a 
system of size Lx L is investigaticd. The results presented 
in Figs. [lj[4]are obtained with L = IOttxo using periodic 
boundary conditions. We expand the wave function as 



*(x)=]Te^ 
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(3) 



and truncate the expansion, \n x \, \n y \ < TV; empirically 
we found that N w 100 is large enough for our purposes. 
For the vector potential we use the gauge A = (0, Bx, 0) 
with a discontinuity at X = ±L/2. The quantity to be 
calculated is the (spin dependent) local density of states 
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FIG. 1: Local density of states (arbitrary units) as a function 
of energy for the Hamiltonian |2| in the absence of a poten- 
tial V(x.) but in the presence of a magnetic field B = Bq. 
The energy is measured in units of Eq — fep/io- The inset 
shows the energy E as a, function of the Landau level index n 
(counting starts at E — 0). 
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FIG. 2: Local density of states in the presence of a Gaussian 
potential with Vo = —2Eq and a = 3a;o- The curves corre- 
spond to magnetic field strengths B/Bo — 0,0. 2... 2 (from 
bottom to top). With the parameters given in the text the 
energy scale is Eg « 33 meV and the magnetic field scale is 
Bo ~ 6.6 T. The potential has a depth of 66 meV and a width 
of 30 nm. 



defined as 
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The spinors 'FaCs, x) are eigenfunctions and e\ are eigen- 
values of the Hamiltonian (pi) . To evaluate the density of 
states at x = we apply a method based on an expansion 
in Chebyshev polynomials following Ref. [23 . Instead of 
a sequence of (^-functions the Chebyshev expansion pro- 
duces a broadened density of states, the broadening be- 
ing determined by the number of polynomials kept in the 
expansion. 

The accuracy of the method is demonstrated in Fig. [T] 
where we show the local density of states in the presence 
of a magnetic field (B = Bo) but in the absence of the 
potential V(x). At low energies one observes a sequence 
of Landau levels at positions E n = ±.Eo\j2nBjBo as ex- 
pected from the analytical calculations. The spectrum is 
reproduced with a very high accuracy as it is shown in 
the inset of Fig. [T] Due to the truncation of the Cheby- 
shev expansion the Landau levels in our figure have a 
finite width. As a consequence only a few discrete peaks 
are seen in the low energy region of the density of states, 
whereas at larger energy the Landau levels overlap such 
that one observes the linear density of states of the Dirac 
Hamiltonian. 

Figure [2] depicts the local density of states in the 
presence of a Gaussian potential with Vq = —2Eq and 
a = 3a;o ■ The magnetic field varies from B = to 
B = 2B in steps of AB = Q.2B . For clarity we use 
offsets for the different magnetic field strengths. With 
xq = 10 nm and vf — 5 x 10 5 m/s the energy and mag- 
netic field scales are Eq rj 33 meV and Bo w 6.6 T. With 
our choice of parameters we are close to the values given 
in Refs. [T7] and [TH]. For zero magnetic field (lowest 



curve), the potential shifts the minimum of the density 
of states to a lower energy, and the density of states is 
no longer symmetric around the minimum. In the pres- 
ence of a magnetic field one observes several Landau lev- 
els. However, the negative branch of the Landau levels 
is suppressed. Only for large magnetic field peaks that 
correspond to negative Landau level index n reappear. 

In order to obtain a qualitative understanding of these 
findings we present now a semiclassical analysis of the 
effect. We start with the classical equations of motion 
for wave packets formed by the eigenstates of the Hamil- 
tonian ([2|, 



v = ±upk 
hk = -ev x B- W(x); 



(6) 
(7) 



here v is the velocity, k is the direction of k and the plus 
and minus sign correspond to particles and holes respec- 
tively. In the absence of the potential V(x) the magnetic 
force constrains the trajectories to circles with cyclotron 
radius r c — hk/eB. Using the Bohr-Sommerfeld quanti- 
zation condition 



I = w p • dx = 27rh(n + 7), p = Kk — eA 



(8) 



with 7 = the correct Landau level spectrum is recov- 
ered. Generally, 7 is the sum of a Berry phase and a 
Maslov contribution which cancel each other in the case 
of massless Dirac fermions, see for example Ref. [21] for 
a recent discussion. 

In the presence of the electrostatic potential V(x) there 
is a competition between magnetic and electric forces. 
Considering for example a circular motion around the 
origin we find a modified cyclotron radius, 



hvpk 



ev F B±V'(r c )' 



(9) 
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FIG. 3: Semiclassical analysis of the Landau levels for B = 
Bo in the presence of a potential V( x ) • The upper curve shows 
the action I/h for closed trajectories starting and ending x = 
(right axis), the lower curve shows the local density of states 
AT(E,x = 0) (left axis). The depth and width of the potential 
V(x) are the same as in Fig. [2] 
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FIG. 4: Spin resolved density of states Af s (E, x = 0) for the 
same parameters as in Fig. [3] 



For a Gaussian potential with negative Vq the deriva- 
tive V'(r c ) is positive and thus the potential reduces the 
cyclotron radius for electrons (due to the extra force di- 
recting towards the center of the cyclotron orbit) but 
enlarges the cyclotron radius for holes. An analogous be- 
havior is found for the action integral I: for a given k, 
the potential increases the action of closed loops for holes 
and decreases the action for the electrons. 



The upper part of Fig. [3] depicts numerical results for 
the action / (divided by K) for trajectories that start and 
end at x = as a function of energy. The trajectories are 
closed, in general they are not periodic but form rosette- 
like orbits. Clearly for electrons (E > Vo = —2Eq) the 
action grows much slower as a function of energy than for 
the holes just as we argued before for circular trajecto- 
ries. In the lower part of Fig. [3] we show again the local 
density of states at x = 0, but compared to Fig. |2]we in- 
creased the number of Chebyshev polynomials such that 
the peaks become sharper. Due to this improved resolu- 
tion the density of states below E rj —2Eq is no longer- 
smooth but also shows a pronounced peak structure. Fur- 
thermore, with exception of the two peaks at E s=s — 2E$ 
and E « — 3E , the Bohr-Sommerfeld quantization con- 
dition accurately reproduces the peak positions. 

Figure [4] shows the spin resolved density of states for 
the parameters of Fig. [3] The central peak is fully spin- 
polarized which is not a surprise since already in the ab- 
sence of V(x) the n — Landau level is spin-polarized. 
More surprisingly the electrostatic potential V(x) splits 
the higher Landau levels into two spin polarized peaks, 
an effect for which we do not have a semiclassical expla- 
nation at the moment. 

In summary we investigated the Landau levels in the 
surface states of a topological insulator. In order to 
mimic the influence of an STM tip on the local density 
of states we included an electrostatic potential in the 
description. For a realistic set of parameters the local 
density of states behaves very similar to what has been 
observed experimentally, namely the negative branch of 
Landau levels appears to be suppressed. The origin of 
the effect is the widening of the cyclotron orbits of the 
holes due to the electric force. We notice that similar 
STM studies exist also for the Landau levels in graphene, 
where the low energy physics is governed by Dirac cones 
as well. In a study of graphene on graphite, where metal- 
lic graphite can screen the electrostatic effects due to 
the STM tip, both branches of Landau levels have been 
observed^. On the other hand in graphene on insulat- 
ing SiC>2 substrates there is no such screening and only 
one branch of Landau levels is observed experimentally^^, 
similar to what was found for Bi2Se3. This suggests that 
in both cases the electrostatic field due to the STM tip 
is the origin of the asymmetric STM spectra. 
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